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Abstract 

Euclidean n-component <f)^ theories whose Hamiltonians are 0(n) symmetric except for quadratic 
symmetry breaking boundary terms are studied in the film geometry R ~ x [0, L] . The boundary 

terms imply the Robin boundary conditions dn4>a = 4>a at the boundary planes 23j=i at 2: = 

0(7) 

and !Bj=2 at z = L. Particular attention is paid to the cases in which rrij of the n variables Ca 
associated with plane *8_, take the special value Cm -sp corresponding to critical enhancement 
while the remaining ones are larger and hence subcritically enhanced. Under these conditions, the 
semi-infinite system with boundary plane !8j has a multicritical surface-bulk point, called rrij- 
special, at which an 0{mj) symmetric critical surface phase coexists with the 0(n) symmetric 
bulk phase, provided d is sufficiently large. The L-dependent part of the reduced free energy 
per cross-section area behaves asymptotically as Ac/i'^~^ as L ^ 00 at the bulk critical point. 
The Casimir amplitudes Ac are determined for small e = 4 — d in the general case where mc,c 
components 4>a are critically enhanced at both boundary planes, mc,D + "Id,c components are 
enhanced at one plane but satisfy asymptotic Dirichlet boundary conditions at the respective 
other, and the remaining tmd.d components satisfy asymptotic Dirichlet boundary conditions 
at both 25 j. Whenever mc.c > 0, the corresponding small-e expansions involve, besides integer 
powers of e, also fractional powers e*""^^ with fc > 3 modulo powers of logarithms. Results to 
order e^^"^ are given for general values of mc,c, ?tic,d + JTId.c, and mD,D. These are used to 
estimate the Casimir amplitudes Ac of the three-dimensional Heisenberg systems with surface 
spin anisotropies for the cases with (mc,c, "1c,d + ?tid,c) = (1, 0), (0, 1), and (1, 1). 
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1. Introduction 

The theory of finite size effects at continuous phase transitions [1] has advanced con- 
siderably during the past decades [2-6]. ^ There is growing evidence that the finite-size 
scaUng functions it involves become independent of microscopic details provided appro- 
priate variables (scaling fields) are chosen and certain nonuniversal metrical coefficients 
are fixed [7,8].^ Universal finite-size properties generally depend on (i) all gross bulk 
properties that determine the bulk universality class (such as space dimension rf, type of 
order parameter, and gross features of the interactions), (ii) gross features of the bound- 
ary (such as large-scale boundary conditions and other properties that determine their 
boundary critical behavior), and (iii) geometric properties such as shape or curvature. 

Here we shall be concerned with d-dimensional systems in a oo'^~^ x L film geometry. We 
have in mind systems of classical n-vector spins whose bulk critical behavior is described 
by a standard 0{n) invariant (f)^ model. Generalizations of this model to films have been 
studied in a number of papers [9-12]. A finite-size quantity of primary interest is the free 
energy density. The total free energy per cross-section area A and thermal energy ksT 
of such systems can be decomposed as 

fL{T) = hm = L MT) + /,(T) + /,es(i; T) (1.1) 

j4^oo AKbJ- 

in the limit A — > oo, where /& and fs are the reduced bulk free-energy density and the 
reduced surface excess free-energy density, respectively, both of which arc independent of 
L. The remaining X-dependent contribution is the residual free-energy density fies{L; T). 
We have kept only the thermodynamic field T, assuming the absence of magnetic fields. 

For films whose interactions are homogeneous throughout the sample except for local 
modifications near the top and bottom planes of the film, depends merely on the 
bulk interaction constants. By contrast, fs and f^^s additionally depend on the modified 
( "surface" ) interaction constants at both boundary planes. According to finite-size scaling 
theory, the singular part of f,-cs{L;T) should have a scaling form which directly at the 
bulk critical point T = T^^oo, yields the asymptotic large-L behavior 

/,ts(L;r,,eo)«A^^^i~(''-'^ . (1.2) 

The associated amplitudes A[f\ called critical Casimir amplitudes, depend on the bound- 
ary conditions p that hold asymptotically on large length scales. They have been com- 
puted by means of the mean-field approximation, approximate renormalization-group 
(RG) methods, exact model solutions, computer simulations, and conformal field theory 
for periodic (p = P) and antiperiodic (p = AP) boundary conditions, as well as a variety 
of large-scale boundary conditions one encounters for lattice systems with free boundary 
conditions [9,10,13-20]. 

The latter large-scale boundary conditions may be characterized via the universality 
classes pertaining to the boundary critical behavior at the bulk critical point of semi- 
infinite systems bounded by a pair of surface planes of the corresponding kinds at 2; = 
and z = oo.li the Hamiltonian is 0{n) symmetric, three distinct familiar types of such 
"surface universality classes" may be distinguished: those associated with the ordinary, 

^ Extensive lists of references can be found in the review articles [2—6] . 
For recent discussions of these issues, see references [7] and [8] . 
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special, and extraordinary transitions [21-23]. Which one of these transitions occurs, 
depends on whether the enhancement of the pair interactions at the surface relative to 
those inside the film is weaker than, equal to, or stronger than a threshold (the "critical 
enhancement"). For those values of n and d for which all three kinds of surface transi- 
tions are possible in semi-infinite geometry, we can specify the boundary conditions p 
by a pair (a, b) of surface universality classes a,b = ord, sp, and ex. Since for supercrit- 
ical enhancement, a transition to a surface ordered, bulk disordered phase occurs at a 
temperature Ts^oo > Tc.oo, the symmetry is spontaneously broken at the extraordinary 
transition. Therefore, boundary conditions (a, b) with a = ex or 6 = ex differ qualitatively 
from those involving only a,b = ord and sp inasmuch as the former break the internal 
symmetry 0{n) down to 0{n — 1), whereas the latter preserve it. Symmetry breaking 
large-scale boundary conditions are also encountered for arbitrary (sub- or supercritical) 
enhancement when boundary contributions linear in <p are present in the Hamiltonian. 
Such a situation is gencrically encountered in the scalar case of one-component and bi- 
nary mixed fluids bounded by walls, i.e., the so-called normal surface transition [22-25]. 
The explicit symmetry breaking (by linear boundary terms in the Hamiltonian) one has 
there is physically distinct from the spontaneous symmetry breaking that occurs in the 
surface ordered, bulk disordered phase and hence is present at the extraordinary transi- 
tion. Nevertheless, the asymptotic boundary conditions are the same in both cases; they 
correspond to fixed nonzero values of the order parameter at the surface. 

In this paper we will be concerned with the Casimir effect under symmetry breaking 
boundary conditions of a different kind. Our aim is to generalize previous work on the 
Casimir effect in models with 0{n) symmetric Hamiltonian [10] by including quadratic 
symmetry-breaking boundary terms. Our motivation is the following. The spontaneous 
breaking of the continuous 0(n) symmetry that the extraordinary transition involves re- 
quires the surface dimension d — 1 to be larger than 2. This means that at d < 3 neither 
the extraordinary nor the special surface transitions are possible for 0{n) symmetric 
Hamiltonians and free boundary conditions. ^ By contrast, normal transitions are possi- 
ble whenever d exceeds the lower critical dimension d* = 2 of the 0{n) bulk model. Thus, 
for three-dimensional isotropic n-vector models the possibility of a special transition at 
finite critical enhancement does not exist. ^ However, in systems whose interactions are 
known to be isotropic in the bulk (such as the isotropic Heisenberg model) it is quite 
common that surface spin anisotropics occur, which entail quadratic symmetry-breaking 
boundary terms in the Hamiltonian. When they are allowed, the picture changes and 
becomes much richer. Suppose the boundary terms at plane are isotropic in an mj- 
dimensional subspace with 1 < rrij < n. Then a surface transition breaking this 0{mj) 
symmetry can occur if the surface coupling associated with this easy manifold is suffi- 
ciently strong and d—1 exceeds the respective lower critical (bulk) dimensions 2 and 1 for 
continuous (m^ > 1) and Z2 symmetry, respectively. Consequently, anisotropic analogs of 



^ The case n = 2 is special since a surface phase with quasi-long-range order is possible at d = 3. 
Thus a surface-bulk multicritical point at which a line of surface transitions of Kosterlitz-Thouless type 
reaches the bulk critical temperature may be anticipated. We will not embark on an investigation of this 

possibility. 

In a three-dimensional semi-infinite lattice model whose spins are coupled by ferromagnetic nearest- 
neighbor bonds ^1 of strengths K\ and K depending on whether both sites x and x' belong to the 
surface layer or not, the entire line of extraordinary transitions including the surface-bulk multicritical 
point associated with the special transition should degenerate into a point with {Ki/K)sp = 00. 
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the special transition, so-called -special transitions [26,27], occur in the corresponding 
semi-infinite systems when the enhancement of this coupling takes a special value while 
the remaining surface couplings are weaker. For d = 3, all cases of critically enhanced 
easy-axis spin anisotropies — namely, those where either mi or m2 or else both are unity 

— arc of interest. 

One could also consider the cases where some or all remaining surface couplings are 
stronger (i.e., supcrcritically enhanced). Following Rcf. [27], we shall exclude this possi- 
bility for the sake of simplicity in all our explicit calculations below. However, we shall 
investigate the general case of a (d = 4 — e)-dimensional film whose surface interactions 
are such that the respective semi-infinite systems with surface planes QSi and 052 would 
undergo mi- and m,2-special surface transitions of this kind with < m,j < n at Tc,oc- 
The weak (not critically enhanced) surface interaction constants are driven to a fixed 
point at which those c;oniponents of with which they are associated satisfy Dirichlct 
boundary conditions. Extending our previous analyses [28,29], we shall determine the 
series expansions of the corresponding Casimir amplitudes A^f ■* in e = 4 — d to order 
e3/2_ Whenever components exists whose surface interaction constants are critically 
enhanced on both boundary planes, these series expansions involve fractional powers e*^/^ 
with A: > 3 (modulo powers of logarithms) besides integer powers of e. 

The remainder of this paper is organized as follows. In Section 2 the model is intro- 
duced. Section 3 sets the stage for the subsequent RG analysis, providing some necessary 
background on the model's renormalization, on anisotropic special transitions in semi- 
infinite systems, their fixed points, and the modified RG-improved perturbation theory 
required in those cases where zero modes arc present at Tc,oc for finite L in Landau 
theory. Section 4 describes our calculation of critical Casimir amplitudes and presents 
our results. Section 5 contains a discussion of our main findings and concluding remarks. 
Finally, there are two appendices with technical details. 



We consider a 0** model for an n-component order-parameter field 0(a;) = {(t)a{x)] a = 
l,...,n) defined on the d-dimensional slab 03 = W'-~^ x [0, L]. Introducing Cartesian 
coordinates Xj, we write the position vector x = {x\, . . . ,Xd) as x = {y, z), where y = 
(xi , . . . , Xd-i) is the (d— l)-dimcnsional lateral component, while z = is the coordinate 
across the slab. It is understood that periodic boundary conditions are chosen along all y 
directions, so that the boundary 903 consists of the union !Bi U!B2 of the two hyperplanes 
z = and z — L. 

We assume that long-range interactions are either absent or can be ignored, and that 
the modifications of the interactions in the vicinity of the boundary planes 05 1 and 052 
are short ranged. Under these assumptions a mesoscopic description in terms of a local 
field theory with a Hamiltonian of the form 



is possible [22,23]. Here the volume and surface elements are given by dV = d'^'x and 
d^ = d'^^^y, respectively. We orient the boundary 903 such that the normal n points 
into the interior of 03. 



2. Model 




(2.1) 
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The bulk and surface densities £(u and £.a<xs depend on and its spatial derivatives. 
We choose for the former the 0(r?) invariant expression 

1 . " 

= 2 E(^'^«)' + i + 5 ' (2-2) 

a=l 

where (j) = \(f)\ is the absolute value. The surface density is taken to consist of general 
terms quadratic in that break the 0{n) symmetry of £sjj; we write it as 

1 " 

i^dm[4>] = o XI ca(a;)0a(a;)^ , (2.3) 

where the surface enhancement variables Ca{x) are allowed to take different, yet position- 
independent, values on the planes *Bi and *B2. In other words, we assume that 

Ca{x) = c^^ forcceOSj. (2.4) 

From the boundary terms of the classical equations of motion (implied by 671 = 0), 
one obtains the boundary conditions 

a„0„(X) = c«^aW forXeOS,-. (2.5) 
These mcsoscopic (Robin) boundary conditions hold beyond Landau theory in an op- 
erator sense — i.e., inside of averages — for the bare regularized theory [22,23]. In our 
actual calculation described in Section 4.1 we shall use dimensional regularization, al- 
though we shall occasionally comment on how the required counterterms depend on a 
large-momentum cutoff if such a regularization scheme were used instead. 

Prom the boundary conditions (2.5) and the surface density (2.3) it is clear that (f)^ 
will get strongly suppressed at 5Bj when Ca takes large positive values. In the limit 
clP Qo, Eq. (2.5) turns into a Dirichlct boundary condition for on *Bj. These 

observations tie in with the known fact that the deviation of —Ca from a reference 
value is a measure of how much the pair interaction on is enhanced relative to its 
bulk counterpart [21,22]. For our subsequent analysis it is also important to remember 
that the mesoscopic boundary condition (2.5) must be carefully distinguished from the 
large-scale boundary conditions mentioned in the Introduction. Just as other coupling 
constants of the mesoscopic theory, the variables Ca^ depend on the length scale to which 
one has coarse-grained. Integrating out the degrees of freedom up to a larger length scale 
generally produces modified values of the interaction constants. In a field-theoretic RG 
approach, their renormalized analogs become running variables. The asymptotic large- 
scale boundary conditions are associated with RG fixed points [22,23]. To work out the 
details for our generalized model with surface spin anisotropies, we first must recall some 
background on semi-infinite systems involving quadratic symmetry breaking boundary 
terms. 

3. Background 

3.1. Free propagators and boundary conditions 

Consider the {N + Af)-point cumulants involving N interior points Xj ^ and M 
boundary points £ 52J in d = 4 — e dimensions, 



5 



/ ^ \ cum 

GiSM(=^l'---'^M) = {]l'^aA^j)]l<l^0dXk)) . (3.1) 

\ j=l fe=l / 

In the disordered phase, the Fourier y-transform 

Gl,„(p;2i,^2) = J d''-^yGLAyi,zi;y2,Z2)e-'P-y'' (3.2) 

of the a = (3 element of the free propagator GL,a{xi;x2) Safj satisfies the differential 
equation 

{-dl+p''+T)GL,a{p;ZuZ2)=5{z^-Z2), (3.3) 

subject to the boundary conditions 

-d,GLAP;z,Z2)\,^L = C^a^GLAP;L,Z2), (3.4) 

where Xj = {yj,Zj) and yi2 = Vi — y2- The solution can be expressed in terms of 

the eigenvalues [fcr(ci,C2)]^ and cigcnfunctions (fr{C,\ci,C2) of the operator —9^ on the 
interval [0, 1] obeying the boundary conditions 

(/?^(0|ci, C2) = ci (/?r(0|ci, C2) , 

-i^^(l|Ci, C2) = C2 ¥'r(l|Cl, C2), (3.5) 

where ci,C2 € [0, 00) and C,= z/L [30,31]. This gives 

iPr{zilL\Lc^^\ Lc'i'^) ^;{z2/L\Lc^^\Lc^^^) 



(3.6) 



The eigenfunctions ipriQ = z/L\c\, C2) are phase shifted cosine functions oc cos(krQ-\-'&r) 
whose phase shift iJr = "d-r (01,02) follows from the boundary conditions at z = 0. The 
boundary condition at z = L yields an equation for the eigenvalues = (kr/L)^, which 
for general values of ci and C2 is transcendental. Details can be found in Refs. [30] and 
[31], but will not be repeated here since we shall need the explicit expressions for the 
free propagators GL,a(p; zi, Z2) and the eigenfunctions only for the special choices o!i^ ~ 
0,00. For those, the free propagators satisfy Dirichlet (D) or Neumann (N) boundary 
conditions on *Bj, depending on whether Ca^ = 00 or Ca^ = 0. The eigenfunctions and 
eigenvalues then simply become 





= (^(C|o,o) = 


= ^J2-5rfi cos(fcrC)) 


kf 


= rw, 




= (p(c|oo,oo; 


) = \/2sin(fc^C), 


kf 


= rw, 




= <^(C|0,oo) 


= •\/2cOs(fcrC)) 


kf 


2r+ 1 
2 


^P'^HC) 


= V5(C|oo,0) 


= •\/2sin(fcrC)) 


kf 


2r+ 1 
2 



r = 0, l,...,oo, 

r e z, 

r = 1,2,. ..,00, 

r = l,2,...oo. (3.7) 
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The respective free propagators can be written in terms of G^'{p;z), the Fourier 
y-transform of the free bulk propagator 

GS^(^I^) = / 4-7 = i'^^r'^' {f/xY^'^^ K^,-2y2{xVf) (3.8) 

J q y \ I 

in d dimensions, where J^^^ = (27r)~^ J^^ ^^Q is a convenient short-hand. One finds 

oo 

4^'^)(p; zi, Z2) = ^(-l)^d-^AB) _ ^2 + 2jL) 

- (- 1)**n5bn + ^2 + 2jL)] , A,B = N,D. 

(3.9) 

These results can be obtained directly by means of the method of images. Alternatively, 
one can start from Eq. (3.6), substitute the eigenfunctions (3.7), and then use Poisson's 
summation formula (see, e.g., Eq. (4.8.28) of Ref. [32]) 

1 °° ^ r°° Ah 

1 Y: f{rn/L) = Y: / -me^'^'' ■ (3.10) 

r=-oo j=-oo''-°° 

Closed- form expressions for G^'^^ (p; zi,Z2) will not be needed but may be gleaned from 
Eq. (B5) of Ref. [8]; its right-hand side, with the frequency variable w appearing there set 
to zero, gives GL,a{P', zi, Z2\T,Ca\ca^) for general non-negative values of Ca^ and Ca^. 



3.2. Renormalization by expansion about the isotropic special point 



"(1) 

Upon expanding in the variables cd , we can work with the free Neumann- Neumann 
propagator G^^'^\xi;x2)- This has integrable ultraviolet (UV) singularities at coin- 
ciding points originating from its bulk contribution G^^a^i — X2) and additional ones 
localized on the boundary planes Q3i and *B2 due to two other summands in Eq. (3.9). 
As is well known, these UV singularities induce divergences in the Fcynman integrals for 
^a^'^^M ^hich can be absorbed in a systematic fashion by appropriate reparametriza- 
tions when d < 4. In addition to bulk countcrtcrms, countcrterms with support on the 
boundary planes are needed [22,23,33-37]. They can be chosen to correspond to the 
reparametrizations known from the bulk and semi-infinite analogs of the model. We use 
the bulk reparametrizations 

= Zy^0ij, St = T - Tc,oo = Zr IJ-'^T, U Nd = jJ." u, (3.11) 

to introduce a renormalized field 0/?, renormalized variable r, and coupling constant u, 
where ^ is an arbitrary momentum scale. Following Ref. [29] and [8], we choose the factor 
that is absorbed in the renormalized coupling constant as 



Nd= 2r(3-d/2) _ 1 



l + l^^^±Ml^e + 0(e^) 



(3.12) 



(d-2)(47r)'^/2 167r2 

Owing to the presence of the 0(n)-symmetry brciaking boundary terms, the usual sur- 
face reparametrizations [37,22] must be generalized as expounded in Refs. [27] and [38]. 
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Accordingly, we introduce renormalized surface enhancement variables Ca ' and renor- 
malized boundary fields (plgrg via 

n 

Scfi^ =liY,Za0cf, Sc^^ = c^^-c,p, ^\a<B = {Z^Z.f/^ ct>\§^ , (3.13) 

13=1 

where ^\dvs means <p, taken at a boundary point X, and Zap has the form 

Zai3 = [Zc + (n<5a/3 - l)Zf]/n. (3.14) 

We fix the rcnormalization factors Z^, Zt, Zu, Z\^ Zc, and Zj by minimal subtraction of 
poles in e = 4 — d. Then they agree with the results to order given in Eqs. (3.142a)- 
(3.142c), (3.66a), (3.166b) of Ref. [22] and Eq. (49) of Ref. [27]. 

It is convenient to transform to enhancement variables that transform orthogonally 
under the RG. Let U = (Uafj) be an orthogonal matrix with = n^^/'^. Upon intro- 
ducing the symmetric deviations Sca^ and anisotropy variables fa^ via 

r/(j) fora = l,...,n-l, 

n 

V C/„fl ^cL^'^ = < _ " , (3.15) 
^1 ^ ^ |5c« = ^4^Vn fora = n, ^ ' 

I 0=1 

as well as anisotropy operators 

n 

Aa = J2Ua0^}, a = l,2,...,n-l, (3.16) 

0=1 

(which we do not further specify; for explicit expressions, see Ref. [27]), the boundary 
operators describing the deviations from the isotropic special point become 

a=l 0=1 

Further, the corresponding linear combinations of the renormalized enhancement vari- 
ables Ca ■* are related to these bare variables in the simple multiplicative fashion 

5c^^^=liZj^\ f^J^ =liZff^j\ a^\,...,n-\. (3.18) 

It follows that the rcnormahzcd fimctions G^^'^^' = ^-(^+^)/2^-A^/2G(iv,M) ^^^iose 
component indices ai , . . . , (3m we suppress) satisfy the RG equations 



' M+N M ' 
T^u- + — 2 — ^'^ + Y ^1 



G^f '^^ = (3.19) 



with 

2 2 n-1 

= nd^ + Mu - (2 + Tlr)dr - (1 + rie) E c^'^^cO) " (1 + ^/) E E ■^^'^^/<^> ' (3-20) 

3=1 j=l a=l 

where the beta function /?„ and the exponent functions riw are defined through the fi- 
derivatives 9^|o at fixed bare interaction constants, 

/3u = m9^|ow= -[£ + %(«)]«, Vw{u) = iJ-d^lo^nZyj, w = <j),T,u,l,c,f. (3.21) 



The characteristics of these RG equations define us running interaction variables u{£), 
f(£), c(£), and fa\t) describing their flow under changes ji ii£. The flow equations 
of u{£), f{i), and c{£) and the respective initial conditions can be found in Ref. [22, 
Eq. (3.79a-c)]; those of the fi^\i) read 

e^/i'\e) = -[^ + Vfiu)]fi'Hi), fi'H^)=fL'^- (3.22) 

The asymptotic behaviors of the isotropic enhancement c{£) and the anisotropics fa^ {£) 
in the large-length-scale limit ^ — > are governed by the familiar (surface-bulk) crossover 
exponent 

$ = + 77e(w*)] (3.23) 
and the anisotropy crossover exponent (cf. Ref. [27]) 

^ = iy[l + r]f{u*)], (3.24) 

where u* = 0{e) is the nontrivial root of /3„ (infrared-stable fixed point). One has c{£) ~ 
£-^/''c and fi^\£) ~ 

As is well known, renormalization of the bulk and surface free energy densities fb and fg 
requires, in addition to the counterterms implied by the bulk and surface reparametriza- 
tions (3.11) and (3.13), also additive counterterms. Owing to the latter, the corresponding 
RG equations are inhomogeneous. However, the additive counterterms can be chosen in a 
way that their contributions drop out in the residual free energy /res [22] , which therefore 
satisfies the homogeneous RG equation 

= 0. (3.25) 

Solving this at the bulk critical point T = Tc^oo (and neglecting corrections to scaling), 
we arrive at the scaling form 

/res(i; Tcoo) « D {c^^^ L"^ ' , c^^^L^/"- {/^^U*/''}) L-^'^-^\ (3.26) 

The restriction of the function £'(ci,C2; {fa^}) to vanishing values of the anisotropy 

scaling variables f^^ = fa ^'L^^^ is precisely the sealing function Z)(ci,C2) determined 
to order e in Ref. [12]. If we additionally set Cj to the fixed-point values and oo of 
the isotropic special and ordinary transitions, we must recover the respective Casimir 

amplitudes A(^P'^p\ ^i^^^ord) ^ ^(ord.ord). 

I?(0,0;{0}) forp = sp,sp, 

D{oo, oo; {0}) for p = ord, ord, (3.27) 
D{0, oo; {0}) for p = sp, ord. 

The expansions of ^(^''^'"'''i) ^ and of a[^'''^''' and a[^^'°'^'*\ were first determined to 0(e) in 
Refs. [34] and [9] , respectively. That the small-e expansion of A^^''"'''''' involves fractional 
powers of e modulo logarithms beyond 0(e) was shown in our work with Shpot [28] and 
Ref. [29]. 

To proceed we must first recall some background related to the anisotropic special 
transition in semi-infinite systems. 
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3.3. Anisotropic special transitions 



In order to characterize m-special transitions in semi-infinite systems, it will be helpful 
to introduce the (zero-field) surface susceptibilities 

/ dAGl^;%{X„X,)=Xjk;a,Sa,^,, (3.28) 

J'Bj 

where Xk is an arbitrary point on *Bfe. The function Xjk-.a describes the response of 
(j)a{Xj) to a surface magnetic field h'^''^ = h^^^Sa localized on and oriented along the 
a- axis. 

At an TOi-spccial transition of a scmi-infinitc system with surface plane *Bi, mi of the 
n enhancement variables Ca ■* of this plane are equal and take a special value Cmi-sp while 
the remaining ones are larger, i.e., are subcritically enhanced relative to this value Cmi-sp- 
Let us label the m\ "easy axes" spanning the subspace in which the orientation of the 
order parameter is energetically favored on Q3i as ai^^ and the remaining ones pertaining 
to the local "hard axes" as oi-^^ ■ In order that an mi-special transition occurs in the 
semi-infinite system with surface 05 1 (in sufficiently high dimensions d), we must have 

c^a' = Cmi-sp for a e {a'^J-^} , 



c 



> c^i-sp for a e {ai^^} . 



(3.29) 



At such a transition, the local surface susceptibilities Xii;ali,=oo ^^^^ temperature 
singularities of the form 

fr'^^'i foraG{aW}, 
Xn;.\,=^-\ (3.30) 
ICi.a+Cs.ai"^"." forae{ai')}. 

Here Ci.q and C2,a are nommiversal constants. The e expansions of the surface exponents 
and ^[^if^ coincide with those given in Eqs. (7) and (8) of Ref. [27] for the critical 
exponents denoted 711, e and 'yn,h there, where the parameters ruh and rUg must be set 
to n — mi and mi, respectively. " Since ji^^ > 0, while 7^"^^ < 0, the surface 53i exhibits 
0(mi) surface criticality at the transition point in the sense that the Xii;q|l=oo diverge 
or have a cusp singularity depending on whether a — ai^' or a = al^^-* . That is, 

fO forae{a(i)}, 
I const < 00 tor a e 

where the subscript mi-sp indicates that t = and the conditions (3.29) are satisfied. 

Turning to the renormalized surface susceptibility Xii;a,ii|/,=oo) we note that the RG 
equations (3.19) yield the scaling form 

Xii;a,KL=oo « T-^ii X„ (c(i)t-*, {/.r"*}) , (3.32) 



^ In the case of 7ii'g and 7^^';,, their analogs for semi-infinite systems with surface 032, the surface 
susceptibilities X22,a obviously take over the roles of xii,a in Eq. (3.30), and one must set me = m2 
and m/j = n — 7712 in the cited e expansions. 
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where = — r]^{u*) — 771 (u*)] is a familiar surface exponent of the isotropic special 
transition whose e expansion to order can be found in Eq. (3.156f) of Ref. [22]. Since 
r ~ i, consistency of the crossover scaling form with the asymptotic behavior (3.29) at the 
mi-sp transition requires that the scaling functions with a G {a'^^} have singularities. 

To locate the mi-sp transition, we set Ca ^ = ci^^ for all a e {ai^^}- Fulfillment of this 
condition guarantees that wc have an 0(mi) symmetry in the subspace of the mi easy- 
axes components ai^"* (provided this symmetry is not spontaneously broken). Thus, all 
scaling functions Xa with a e {o;!^^} become identical. For given values of the anisotropy 
variables {fj^^}, we can determine the critical value Cmi-sp of c^^^ associated with the 
mi-sp transition from the condition that limt^o l/xii,acli=oo vanishes at c'^-* = Cmi-sp 
(while keeping the constraint Ca ■* = ci^'' for all a G {cce})- The shift Cmi-sp ~ Cgp must 
have the scaling form 

- Cp = [A«] H^(/,«/A«, . . . , £\/.fi'') (3.33) 

To determine the temperature singularity of the surface susceptibility at the rrii-sp 

transition from the singularity of its scaling function is a rather cumbersome procedure 

_ (1) 

because the RG scheme used so far does not directly yield the power r ''11.= in expo- 
nentiated form. As is shown in Ref. [27] , the required information can be obtained in an 
easier and more direct fashion by formulating the RG against an anisotropic background 
rather than the 0{n) symmetric theory used above. Since the order parameter compo- 
nents (j)i3 with /3 e {Ph^} are noncritical on the boundary planes Q3j in the sense that the 
respective inverse surface susceptibilities l/xjj,i3\L=oo do not vanish at the transition, we 
can set the associated bare enhancement variables c^^ with /? e {Ph^} to the fixed-point 
values c*j.jj = 00 and expand about the mj-sp multicritical points of the corresponding 
semi- infinite systems with surfaces *8j . 

To this end we form, by analogy with Eq. (3.15), the symmetric combination 

^i'^ = ^ T. 4^'^ (3.34) 



m 



and rrij^l orthogonal anisotropy variables k = 1, . . . , nij — l, from the coefficients Ca^ 

with a e {ai"''}. The former couples to the 0{mj) symmetric operator J2ae{a'^''} '^a- 

anisotropics g^^^ couple to analogs of the anisotropy operators Aa for the respective rrij- 

dimensional subspaces, which we denote as ^ Let Qjf^J ~ ^^^Qk^ operators 

one obtains through the substitutions (j)'^ 4'a,R ~ ^^^^a ^ {od^^}- Following 

Ref. [27], we can introduce renormalized variables ci^^ and g^^ by writing the boundary 
contributions to the action as 



jC^d^[Z^<pR] dA = 



(3.35) 
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Here Zg' and Zg are respectively the renormalization factors Zc and Z^ given in 
Eqs. (70) and(71) of this reference, with rrie = rrij and rrih = n — rrij. The superscripts 
(■'^ here serve to remind us that these functions depend via raj on the number j of the 
boundary plane 5Bj. Just as Cgp, the critical value Cmi-sp vanishes in a perturbative ap- 
proach in 4 — e dimensions if dimensional regularization is used. (In a theory regularized 
by means of a large-momentum cutoff A, both quantities would diverge ~ A in d = 4 
dimensions.) 

The RG eigenvalues governing the flow of the running variables ci\l) ~ l-Ve^ ci'^ and 
9k\^) ~ (-~^^'>^ 9^^ are given by 

y« = ! + #(«*), 

y« = l + »7^^)(«*), (3.36) 



where the exponent functions are defined as 

77(^(u)=MS^|olnZe,g, (3.37) 

by analogy with Eq. (3.21). The e expansions of vye^^ and vy'J^ may be found from those 
given in Eqs. (87) and (88) of Ref. [27] (for the exponents denoted 6 and ^ there) by 
setting the parameter mh to n — rrij . 

It should be clear that the generalized asymptotic boundary conditions with which 
we are concerned — namely, those for which nij order parameter components (j)a, ol 

{a'e^}, are critically enhanced on QSj while the remaining n — rrij ones with a G {o^l^^} are 
subcritically enhanced there — are associated with fixed points located at (ci^-*, ci^-*) = 
(0,0) on the zero-anisotropy hyperplane gj^^ = with r = and u = u*. These fixed 
points are unstable in the ci^^ and gj^^ directions. By contrast, surface spin anisotropics on 
5Bj that are restricted to the respective subspaces of hard axes a e {a^^^} are irrelevant. 
Components 0^ whose index a belongs to the set of hard axes, {ajf''}, satisfy Dirichlct 
boundary conditions at the plane *8 j . The behavior near the plane can be obtained from 
the short-distance expansion (cf. Ref. [36] and [22, p. 190ff]) 

<l)0,R{X+nAz) C{Az)[MX)U, (3-38) 

where the operator on the right-hand side means the renormalizcd surface operator 
[4)p{X)]ii = {Z^Zi^e)~^/'^(j)j3{X) at a point X G *Bj. The behavior of the function 
C(Az) at distances 5z small compared to L and ^oo (but large compared to microscopic 
distances) is governed by the difference of scaling dimensions of the operators on the 
left-hand and right-hand sides of Eq. (3.38). We have 

C{Az) « const|A^|<"e'' with ae{n,mj,d) = {rjlfl - 'n)/2, (3.39) 

where 77^"'^ is the surface critical exponent ?7|| ^ of Ref. [27] for m/j = n — ruj. Using the 
results obtained there, one arrives at the series expansions 
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, n-2mj-2 ^ 6 + 5m,- - 2n ^ 2 ^r/ *^3^ 
ae{mj,n,d) = ^ u* + ^ («, ) + 0[{u*) ] 

_n- 2mj - 2 12 + (4 + 5n)n - 4(2n + l)mj ^ 3^ 
~ 2(8 + n) ^ + 2(8 + n)3 ^ + '^^^ 

(3.40) 

Of interest is the case of a (d = 3)-dimensional Heisenberg magnet with an easy-axis 
spin anisotropy at 5Bj (i.e., d = n = 3, rrij = 1). To obtain a rough estimate of the expo- 
nent iTe(3, 1, 3), we evaluate the O(e^) expression (3.40) at e = 1. This yields cre(3, 1, 3) ~ 
—40/1331 ~ —0.03. It is plausible that a small negative value results. ^ However, because 
of our naive extrapolation procedure, this estimate is not very reliable. Unfortunately, we 
are not aware of accurate Monte Carlo estimates of this exponent. Although a number 
of detailed Monte Carlo investigations of 0{n) models on simple cubic lattices [40,41,39] 
were carried out during the past decade, none of these allowed for 0(n)-symmetry break- 
ing surface spin anisotropics. 



4. Casimir amplitudes 



4.1. Calculation of residual free energy via RG improved perturbation theory 



We are now ready to turn to the calculation of Casimir amplitudes. The general case of 
asymptotic boundary conditions we wish to consider should be clear from the foregoing 
section: those for which the renormalized enhancement and anisotropy variables c^J^ and 
g\^^ of boundary plane take their fixed point values c^J^ = g^^ = while all bare 
"hard-axes" surface enhancement variables ■* with a e are set to the value c*^^ = 

00. The conditions ci''^ = gj^^ = imply that all renormalized easy-axes enhancement 
variables Ca^ with a G {ai^^} vanish. In general, we therefore have: (i) nic^c of components 
(pa that are critically enhanced on both planes {a G {ai'^''}n{ai^''}); (ii) mD,D components 
satisfying Dirichlet boundary conditions at both planes; (iii) m^^D and ^ components 
that are critically enhanced on !Bi and *82 but satisfy Dirichlet boundary conditions 
on 5B2 and QSi, respectively. Components that are critically enhanced at *Bj have an 
asymptotic near-boundary behavior of the form specified by Eqs. (3.38)- (3.40) near this 
plane, irrespective of whether they are critically enhanced at the complementary plane 
or satisfy Dirichlet boundary conditions there. 

To determine the corresponding Casimir amplitudes Ac via RG improved perturbation 
theory in d = 4 — e dimensions, we set all enhancement variables to the above-mentioned 
fixed-point values and compute the residual free energy (1.2) at the bulk critical point. 
Whenever mo,c > 0, there are free propagators GL,a = G^^'^^ among those of Eq. (3.4) 
whose spectral decompositions(3.6) involve the eigenfunction ip'^^K Since the associated 
eigenvalue vanishes, we encounter a zero-mode problem of the kind dealt with in Refs. [28] 
and [29]. Its origin is that Landau theory in those cases erroneously predicts a sharp 



® For the analogous exponent (?7y^ — rj) /2 of the {d = 3)-dimensional semi-infinite Ising model's special 
transition, the results = 2.482 and j;^'^ = 1.636 given in Rcf. [39] for the RG eigenvalues y/j and 
of the bulk and surface magnetic fields yield (ryj]^ — »;)/2 = j/^ — 1 — j;^^ ~ —0.154. 
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transition at Tc^oc for films of finite thickness L. The infrared singularities associated with 
these zero modes entail that conventional RG improved perturbation theory becomes ill- 
defined at Tc_oo- As shown in Refs. [28] and [29], this can be remedied by a reorganization 
of RG improved perturbation theory. The strategy is analogous to the one used in the 
theory of finite size effects on continuous phase transitions in systems that are finite in 
all, or in all but one, dimensions [42,43]. Its crux is to split off the zero-mode components 
from the TOc,c fields (f>a{x) with a G {ai^''} H {cki^"*} and integrate out the remaining 
degrees of freedom to construct an effective action for the zero-mode components. To 
this end we introduce the mc, c-component field <p(y) with components 

if^iy) ^ L-^'^ I dzMy.z), a e n {ai^^}, (4.1) 







and decompose 4>[x) as 



, _ /i-i/Va(y)+^a(y,^) forae{aW}n{a(2)}, 
" z) for a i n {af)}. ^^''^ 

Obviously, the field if contains the zero-mode contributions to (j). By orthogonality to 
the r 7^ eigenfunctions lySr^'^"*, we therefore have 



( dz^paiy, z) = for aU a £ {a^}''} n {a^^)}. 
Jo 



(4.3) 



Let us introduce the free energy associated with the field ip, F^, and the corresponding 
reduced area density by 

-^ = -lnTr^e-^[^l, /^(L) = fim (4.4) 

where Tr^(.) means the functional integral J I?'0(.). When mc,c = 0, no zero modes occur 
in the spectral decomposition of the free propagators at Tc^cxd! then F^p and ftp{L) coincide 
with the total free energy F and the reduced area density /l, respectively. Whenever 
^c,c > 0, the effective Hamiltonian, defined by 

is nonzero and yields the additional contribution 

- i^-oo - - i^-oo J lnTr,e---[-l (4.6) 

to the reduced area density 

fL = ML) + U{L). (4.7) 
To compute f^{L)^ we can use a standard loop expansion 

-U{L)A = ( ) + X + 0(3-loops). (4.8) 
^ ' 1 1 

Here the dashed blue lines in the second graph represent free ■(/'-propagators G^^L,a = 
[(3o(-A(^°'^°) + f)Sai3Qo]~^, where A(^°'^°' with Aa,Ba = N,D means the Lapla- 
cian, subject to the boundary conditions (Aq,Bq,) applying to this value of a. Further, 
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the operator Qq projects onto the subspace orthogonal to the ko = modes. In the 
pz representation we have 



+ r) ' 



(4.9) 



as foUows from Eqs. (3.6) and (3.7). 
The one-loop contribution 

1. 



( ; ^ -Trln[Qo(-A(A-B„) ^ f)SM 



(4.10) 



splits into a sum of contributions which may be gleaned from the results for Casimir 
amplitudes of the Gaussian theory under the boundary conditions (A^, Bq,), with A^, Bjj = 
D, N, given in Refs. [34], [9], and [10]. The resulting one-loop contribution to the residual 
free energy density /^,rcs|T=Te oo becomes 



T(d/2)C(d) r , , rl .1 , . 

' 2djd/2i)~l ["*D.D + me,c + {2^-" - 1)(™c,D + m^.c)] , (4.11) 



where (.(d) is the Riemann zeta function. 

The two-loop contribution to f^iL) can be expressed in terms of the integrals 



(i;T) 



dz 



■ GL,-4,,aix; x) Gl,4,,p{x; x) 



as 



(4.12) 



(4.13) 
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In previous calculations [9,10,28,29], only the special cases of the integrals (4.12) with 
a — (3 were encountered. We need these integrals at f = also for a ^ fi. Their 
e expansion to first order is worked out in Appendix A. Upon inserting the results into 
Eq. (4.13) and expressing u in terms of the renormalized coupling constant u, we obtain 



/i/;,res(-^) 



L 



-(d-l) 7tM^ + 0{u)] 



4608 



+ "^D,D) (2 + 3mc,c + 3mD,D) 



-I- - (n -I- 2 - QiTic^c - 6mD,D) + 0(e) 



(4.14) 



We turn next to the calculation of the zero- mode contribution f^^rcs- The one-loop 

f \ 

graph yields to Ticff [</?] a contribution quadratic in (p which changes the coefficient 

of from f to 



dz 



2G(^-N)(^;^) + EGi 



^ ^ + I [i^c,c + 2) /(^'^) (L; f ) + ™d,d (i; r) 

+ (tIcD + "^Dx) 7}'^'°-' (L; f )] . 



(4.15) 
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The required integrals 

(•L 



l[^'''\L;f) = ^Gt'^\x;x), A, B = N, D, (4.16) 
with f = are computed in Appendix A. The results are 



(A,B).. 2i-V'^-3)/^C(3~d) 1 forA = B, 

' ^ ' ' r[(d- l)/2)]cos(^d/2)L'i-2 |(23-'*-l) forA^B. 

(4.17) 



Substituting them into Eq. (4.15), expressing u in terms of the renormalized coupling 
constant u, and expanding in e then gives 

2 

U TT 

tl\t=t,,^ =J2^^ + ^(^)]^ P"^^.^ + 3"'D,D -n + 4 + 0(e)]. (4.18) 
The contribution from the zero mode 

A^(L)= -Q - Q +.... (4.19) 

can now be easily inferred from the results of Ref. [29]. To obtain the one-loop contri- 
bution to the residual free energy density f^^resiL) we must simply set n = rric^c in its 
Eq. (4.28) and replace the variable it involves by the quantity given in Eq. (4.18). 
The two-loop term follows in the same manner, except that an additional sign change 
must be made (as can be seen from Eqs. (4.23) and (4.27) of Ref. [29]). Hence we find 

with 

Ad = -{47T)-''/^r{l-d/2), (4.21) 

where the tl on the right-hand side of Eq. (4.20) are to be taken at the bulk critical 
point. Since tl\to oo i^ linear in u, its value at u = u* is of order e. This implies that the 
first term in Eq. (4.20) yields a contribution of order e^^^ to the Casimir amplitude. The 
contribution from the second term is of order e^^*^ and hence negligible to the order of 
our calculation (e^^^)- 



4.2. Results 



Before using the results of the foregoing section to determine the small-e expansion 
of the Casimir amplitudes, we first give the results for the Gaussian {u = 0) critical 
(t = 0) theory in general dimensions d > 1. The corresponding Casimir amplitudes 
A(7 Qauss for the cases that d and ton,n components (/)q, satisfy respectively Dirichlet 
and Neumann boundary conditions on both planes *Bj, while the remaining n — d — 
™n.n = ™d,n + fTT-N,!) components satisfy Dirichlet boundary conditions on one plane 
but Neumann boundary conditions on the complementary one, can be inferred from 
equation (4.11). The reason is that in this Gaussian case, we do not have to worry about 
the zero-mode contributions. Upon considering them for r > [i.e., the first term in 
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equation (4.20)] and then taking the limit f — * 0+, one sees that their contributions to 
fiesiL) vanish at T = Tc^oo- Thus the Gaussian Casimir amplitudes are given by 

Ac, Gauss = -^-^^^^^TT^ ["^D,D + mN.N + (2^"'* - l)(mN,D + mu,^)] . (4.22) 

Returning to the interacting case w ^ 0, we can add the results presented in Eqs. (4.11), 
(4.14), and (4.20) to obtain the renormalized residual free energy fres,R{L, T, u). Its value 
at u = u* = 3e/(n + 8) + O(e^) gives us the Casimir amplitude 

Ac = L'^-' fresML, Tc,oc,U*) . (4.23) 

Its small-e expansion becomes 

Ac = ao + aie + 03/2 e^/^ + 0(^3/2^ (4 24) 

with 

TT^ f 7n\ , , 

ao = -—\mc,c + mB,D - j, (4.25) 



^2 r 



ai 



1536 



TOc,c + mD,D 
+ 



7n\f ^ , , , 281n2 2C'(4) 
-j(7.-l + ln(.) + ^-^ 

16 n In 2 (mc,c + 'tid,d)(2 + 3mc,c + 3mD,D) 



225 



n{n + 2 — Gtoc.c ^ Sttid.d) 



3(n + 8) 



and 



03/2 = 



mc.cTT 



288\/3 V 



/ 3wc 



3mD,D - n + 4\ 
n + 8 y 



3/2 



(4.26) 



(4.27) 



where 7^ = — r'(l) is the Euler-Mascheroni constant. 

These expansions generalize previous results for critical Casimir amplitudes for those 
cases in which the 0{n) symmetry is not broken, neither spontaneously nor explicitly 
nor by boundary conditions, by the inclusion of symmetry breaking quadratic boundary 
terms. 

Note that our results (4.24)-(4.27) depend only on the sum of rric.c + tod,d- Since n — 
^0,0 — TiD.D = TicD+'TiD.c) the expansions of Ac to the given order of f **/^ arc insensitive 
to changes of mc,D and mu,c at fixed n, mc,c, and mD,D that preserve the total number 
TTiui = mc,D + m^D,c of components satisfying "unlike" boundary conditions, namely, 
Dirichlet boundary conditions on one boundary plane and those corresponding to critical 
enhancement on the complementary surface plane. The reason for this insensitivity is that 
the one-loop and two-loop graphs contributing to the given order are traces and local 
expressions in position space, respectively. The full theory does not have this property 
and hence in general will give different results for cases of different mc,D with the same 
values of n, mu,D, and ttIui- The obvious reason is that such cases with different values 
of mc,D are physically distinct. As an illustrative example, one can compare the cases of 
(mc,c, mB,D,'rnc,D, mu,c, n) = (0, 1, 2, 0, 3) and (0, 1,1,1, 3) involving critically enhanced 
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easy-plane interactions on 5Bi and critically enhanced easy-axes interactions on both 

boundary planes, respectively. ^ 

The special cases of mD,D = n (Dirichlet boundary conditions on both planes), mc,D + 
fnu,c = n (all components <j)a satisfy Dirichlet boundary conditions on one boundary 
plane and are critically enhanced at the complementary one), and rric^c = n (all compo- 
nents are critically enhanced at both boundary planes) provide important checks of our 
results (4.24)-(4.27). In the first two cases, we indeed recover the series expansions 



Ac(mD,D = n)/n = - 



1440 2880 
2C'(4) 5n + 2 



2 n + 8 



1 — 7b — ln(47r) 
0{e% 



(4.28) 



and 



Ac(mc,D + mD,c = n)/n ■■ 



In' 



In' 



11520 23040 



1 - 7b ■ 



121n(2) 



ln(7r) 



2C'(4) 5n + 2 
T(4) 



7 n + 8 



+ 0(6^), 



(4.29) 



of Krech and Dietrich [9,10] for the amplitudes Aj)^j)/n and A-o^sp/n. Likewise in the 
third case, the resulting series expansion 



Ac(mc,c = n)/n = - — — + e 



1440 2880 

2 



1 - 7b - ln(47r) + 



2C'(4) , 5 n + 2 



+ 



2 n-l-8 



_e3/2 



TT^ /n-l-2 
72V6\n + 8 



3/2 



(4.30) 



agrees with the one obtained for Agp^sp/n in Rcf. [28]. 

For {d = 3)-dimensional Heisenberg systems (n = 3), the following easy-axes cases are 
of interest: 



(i) TOc.c = 


1, mD,D = 


2, TOc,D = 0, TOD,c = 0, n = 3: 


Ac/3 = 


-0.00685 


- 0.00377 e - 0.00572 e^/^ + o(e 


(ii) mc,c = 


0, mD,D = 


= 2, mc,D = 1, mD,c = 0, n = 3: 


Ac/3 = 


-0.00257 


- 0.00132 e-hO(e^), 


(iii) mc,c = 


= 0, mD,D = 


= 1, nic.D = l,mB,c = 1, n = 3: 


Ac/3 = 


0.00171 + 


0.00230 e+0(e2). 



3/2x 



(4.31) 



(4.32) 



(4.33) 

In Table 1, estimates for Ac/n at d = 3 are given for various choices of mc,c, ?71d,d, 
TOc,D, wd,c, and n. In the columns labeled 0{e^) with ct = 0, 1,3/2 we have listed the 



^ For reasons discussed already in the Introduction, the 2-spccial surface transition one is concerned 
with in the first case is expected to occur only when d > 3. By contrast, the 1-special transitions one is 
dealing with in the second case should be possible in d = 3 dimension. 
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values one obtains from the small-e expansions (4.24) by truncating them at this order 
and evaluating them at e = 1. The estimates shown in the column marked as [0/1] are 
the e = 1 values of the Pade approximants Ac = ao/(l — eai/ag). 

Table 1 

Estimated d = ?> values of Ac/n. In the easy-axes cases (rricci^D Di^^) = Ij'i) a-^d (0,n — l,n), 

the numbers of modes with Dirichlet boundary conditions only on one boundary plane satisfy m^D = 
"^D,c = and mc,D + wiD,c = !> respectively. In the cases (n,mc,c,T'JD,D) = (2,0,0) and (3,0,1) it is 
understood that mc,D = "iD,c = 1- The values in the column labeled MC refer to Monte Carlo results 
reported in Refs. [19], [20], [13], and [14] (which agree with each other). 





n 


mc,c 




0(e°) 


0(e) 


0(e=^/2) 


[0/1] 


MC 




1 


1 





-0.00685 


-0.01166 


-0.02243 


-0.02293 






1 





1 


-0.00685 


-0.01166 




-0.02293 


-0.015 




1 








+0.00600 


+0.01282 




-0.04352 






2 


1 


1 


-0.00685 


-0.01109 


-0.01817 


-0.01793 






2 





2 


-0.00685 


-0.01109 




-0.01793 


-0.015(6) 




2 





1 


-0.00043 


-0.00003 




-0.00022 






2 








+0.00600 


+0.01296 




-0.03711 






3 


1 


2 


-0.00685 


-0.01062 


-0.01634 


-0.01522 






3 





2 


-0.00257 


-0.00389 




-0.00528 






3 





1 


+0.00171 


+0.00401 




-0.00502 




n 


— > oo 


1 


n - 1 


-0.00685 


-0.00595 


-0.00595 


-0.00605 




n 


— > oo 





n- 1 


-0.00685 


-0.00595 




-0.00605 




n 


oo 





n 


-0.00685 


-0.00595 




-0.00605 





In cases with m,c,c = (where no zero modes are present in Landau theory at Tc,oo), 
we expect the quality of the extrapolations to be comparable with the original 0(e) 
estimates of Krech and Dietrich [9,10] for the isotropic n > 1 cases with tod,d = n and 
mc^T) = n, respectively. In the Ising case with n = mD,D = 1 they used a Pade [1/1] 
approximant whose additional parameter was fixed by the requirement that the exact 
value Ac{d = 2) = — 7r/48 ~ —0.65 be reproduced in addition to the e expansion to 
first order. This gave the estimate Ac{n = mD,D = 1, d = 3) ~ —0.015, which is in very 
good agreement with Monte Carlo results. We have refrained from giving this value in 
Table 1 in order to have a better comparison with the n > 1 cases, for which analogous 
exact results that could be used for improved Pade approximants are not available to 
our knowledge. 

Obtaining reliable extrapolations in the cases with mc,c > 0, where fractional powers 
of e (modulo logarithms) appear in the small-e expansions, is an even greater challenge. 
To this end a better understanding of these series would be needed. Therefore, reliable 
information about Casimir amplitudes from alternative sources such as IVIonte Carlo 
simulations, exact model calculations, and other approximate RG calculations would be 
very useful. 

We also included estimates for Ac/n in the large-n limit. This limit is known to bo 
equivalent to an appropriate spherical model with 2-dependent constraints [44]. Its exact 
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solution leads to a self-consistent Euclidean Schrodinger equation [45] whose analytical 
sohition is not known for finite L, not even at Tc. oo- Recent numerical work [46] gave the 
result lim„^tx) '-^cid = 3)/n ~ —0.012(1) for open boundary conditions, corresponding 
to the case mo.D = n. According to previous work [28,29,8], the uncertainties of naive 
extrapolations of the small-e expansions become significant as n increases. That our 
extrapolation gives a value at c? = 3 which is roughly 50% of this numerical large-n result 
is therefore not unexpected. On the other hand, we expect our small-e expansion results 
to be in conformity with the largc-n limits in d = 4 — e dimensions since we were able to 
show this for our previously published analogous results for periodic boundary conditions 
[28,29]. 



5. Summary and concluding remarks 



In this work we studied effective interactions induced by thermal fluctuations in a 

medium that is confined by two parallel {d — l)-dimcnsional boundary planes at its 
d-dimensional bulk critical point. These fluctuation-induced interactions manifest them- 
selves through a contribution to the free energy per fllm area A ^ oo that depends on 
the separation L of the confining planes. At the critical temperature Tc. oo of the L = oo 
(bulk) system, the residual free energy densities fres{L;T) introduced in equation (1.1) 
exhibit the asymptotic power-law decays (1.2) in L. The associated fluctuation-induced 
( "Casimir" ) forces per area, 

J^ciL;T) = -kBT^^^^^^, (5.1) 

therefore become long-rang ed at Tc,oo [6,9,10,47]. They decay as {d - l)A[f^L-'', and 

hence are attractive or repulsive depending on whether the amplitudes A[f ^ are negative 
or positive. 

These forces J^c are the analogs of the familiar Casimir forces in quantum electro- 
dynamics (QED) between two ideally conducting grounded metallic parallel plates at 
separation L caused by vaciuun fluctuations of the electromagnetic field [48,49]. Just as 
the latter forces, their thermal analogs J^c depend on gross features of the medium (such 
as space dimension d, number of components n of the order parameter, etc) and of the 
confining plates (such as boundary conditions). In the QED case, the interaction of the 
electrodynamic field with matter (i.e., the metallic plates) usually is taken into account 
only through the choice of appropriate boundary conditions. Hence the problem reduces 
to the study of free field theories in bounded geometries under given boundary condi- 
tions p. Common physically relevant choices are the combinations (D,D), (N,N), (D,N), 
and (N,D) of Dirichlet (D) and Neumann (N) boundary conditions. Massless Gaussian 
field theories have the simplifying feature that when these boundary conditions are im- 
posed (at the mesoscopic scale at which the continuum description applies), they remain 
valid at larger length scales. The reason is that Dirichlet and Neumann boundary con- 
ditions at either one of the boundary planes correspond to infrared stable and unstable 
fixed points of the corresponding boundary field theories, respectively [22,23]. Choosing 
for each component (pa one of these boundary conditions at both planes fixes the 
Casimir amplitude of the Gaussian theory for given d and n uniquely; the corresponding 
values Ac.Gauss are given in equation (4.22). 
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On the other hand, if instead of a Neumann or Dirichlet boundary condition a Robin 
boundary condition 9„c/)q, = Cq with < ciP < oo is imposed at *Bj, then this 
mesoscopic boundary condition does not remain vahd on larger scales but turns into 
an asymptotic Dirichlet boundary condition in the large-length-scale limit [12,22,23]. 

This means that effective scalc-dcpcndcnt Casimir amplitiidcs, depending on the scaled 

0(7) 

argument cd L, result whenever such a Robin boundary condition is involved. 

The study of critical Casimir forces below the upper critical dimension d* = 4 requires 
going beyond the Landau and Gaussian approximations. Then important changes of the 
scenario just described occur. They result from the fact that the (f)'^ interaction shifts 
the Gaussian fixed point with Cq ^ = (all a), at which a Neumann boundary condition 
applies for </> at *Bj on the level of this free theory, to a nonuniversal (cutoff dependent) 
value Csp [see equation (3.13)]. The associated shifted fixed points are the ones of the 0{n) 
invariant theory called "special" . As already recalled in the Introduction, their physical 
significance is to describe the so-called isotropic special surface transitions [22,23], which 
occur (in sufficiently high space dimensions) for critically enhanced surface interactions. 

Accordingly, Neumann boundary conditions do no longer correspond to RG fixed 
points, neither to stable nor unstable ones. The consequences are threefold: (i) Neumann 
boundary conditions that hold on a mesoscopic scale do not persist on longer length 
scales. The asymptotic boundary conditions on large length scales are determined by the 
RG fixed points to whose basins of attraction the values of the surface enhancement vari- 
ables Ca belong. For simple Ising models with nearest-neighbor interactions whose bonds 
have different strengths in the two boundary layers and elsewhere, Neumann boundary 
conditions are expected to lie in the basins of attraction of the fixed points with Cq^ = 00 
at which asymptotic Dirichlet boundary conditions prevail [12]. (ii) Since the analogs of 
the Gaussian fixed points with Ca = — - the special fixed points — are located in the 
space of bare interaction constants at a nonuniversal value Cgp, the mesoscopic boundary 
condition dn4> = Csp(f> one has for critically enhanced surface interactions is also nonuni- 
versal. (iii) Again, this boundary condition does not persist on longer length scales. On 
large length scales <C i,^oo, a characteristic algebraic near-boundary behavior holds in 
this case of critical enhancement, as predicted by the boundary operator expansion (cf. 
section 3.2). 

The main purpose of this paper was to generalize previous investigations of critical 

Casimir amplitudes for systems with n-component order parameters in a film geometry 
by allowing for 0(n)-symmetry breaking quadratic boundary contributions in the Hamil- 
tonian. Aside from the above-mentioned Gaussian results (4.22), our main findings are 
the small-e scries expansions of the Casimir amplitudes Ac gathered in equations (4.23)- 
(4.27) and the d = 3 estimates presented in Table 1. There are several good reasons for 
the generalization we made. 

First, it is not imcommon that surface spin anisotropics occur in systems whose bulk 
critical behavior belongs to the universality class of the 0{n) model. For this reason and 
from a fundamental point of view, there is therefore interest in a study of fluctuation- 
induced interactions in such systems. Second, even in the absence of such symmetry 
breaking boundary terms, one can consider boundary conditions on the mesoscopic scale 
of the continuum field theory that break the 0(n) symmetry. Common ways of breaking 
the internal Z2 symmetry in the scalar {n = 1) case are choices of ±± and ±^ boundary 
conditions, for which the values of the order parameter at the boundary planes Q3j are 
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fixed to nonvanishing values of same or opposite signs. Such boundary conditions are of 
direct relevance for fluctuation-induced forces that are mediated by binary liquid mixtures 
near their consolute point [16,20,47,50-53]. In the continuous symmetry case n > 1, it is 
natural to consider also twisted boundary conditions for which the local order parameters 
(t>\<Sj at the two boundary planes are aligned along different directions but have the same 
fixed magnitude [2,54,55]. Rather than by such inhomogeneous boundary conditions, one 
can also break the symmetry by homogeneous boundary conditions. Thus one can choose 
different Robin boundary conditions for the n components of the order parameter. This 
is precisely the situation we were concerned with in this paper. Just as in the absence 
of 0(n)-symmetry breaking boundary terms, the boundary conditions that hold on a 
given mesoscopic scale do not necessarily hold on larger scales, let alone asymptotically 
in the large-length-scale limit [12,22,23]. Determination of the boundary conditions in 
the large-length-scalc limit requires the understanding of the respective fixed points of 
the associated boundary field theory and the boundary operator expansions. These issues 
were dealt with in some detail in Sections 3.2 and 3.3. 

An important difference between symmetry breaking by the mentioned inhomogeneous 
(±±, ±^ or twisted) boundary conditions and the latter homogeneous ones should be 
realized. The former, which can be realized by symmetry breaking linear boundary terms 
in the Hamiltonian — i.e., surface fields hj = hjhj whose components hj along certain 
directions specified by the unit vectors hj tend to infinity — imply a nonzero order- 
parameter profile m{z) = {4>{y, z)) both above and below the bulk critical temperature 
Tc,oo even in Landau theory. In the case of quadratic 0(n)-symmetry breaking boundary 
terms (or corresponding homogeneous boundary conditions) , the order parameter profile 
vanishes above T > Tc^oo, unless some enhancement variables are supercritically enhanced 
so that in a temperature region above Tc^oo spontaneous symmetry breaking occurs and 
a nonzero profile m{z) results. When this possibility is ruled out because all surface 
interactions are subcritically or at most critically enhanced, the Casimir force that results 
for temperatures T > Tc^oo is entirely due to fluctuations and vanishes in Landau theory. 
In the case of linear symmetry-breaking surface terms (or corresponding inhomogeneous 
boundary conditions), these forces generically vanish neither above nor below nor at 
Tc^oo- Thus, at and above Tc,oo they are not purely fluctuation induced, a property they 
share with the Casimir force below Tc, oo in the case of homogeneous boundary conditions. 
One may question whether the term "Casimir force" is appropriate when a force is not 
entirely due to fluctuations (as the original QED Casimir force is), but it has become 
customary to use it even in such cases [53]. The critical Casimir interactions we studied 
above are all completely fluctuation induced and hence more compelling analogs of those 
known from QED. This provides a further reason for investigating them. 

A fourth reason is that one does not expect special surface transitions to be possible 
in semi-infinite geometry for systems with 0{n) symmetric Hamiltonians and short- 
ranged interactions. The introduction of quadratic symmetry-breaking boundary terms 
corresponding to easy-axes spin anisotropics opens up the possibility of having large- 
length-scale boundary conditions that are governed by anisotropic special fixed points. 
A fifth reason is that the cases with mc,c > (involving zero modes) provide further 
examples of the breakdown of the e expansion for critical Casimir amplitudes. 

A natural question to ask is how the Casimir interactions investigated in this work 
can be checked by experiments and simulations. Previous experimental work on thermo- 
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dynamic Casimir interactions has focused on fluid systems [50-53,56]. The advantage of 
fluid systems is that the number of degrees of freedom of the medium (fluid) between 
macroscopic objects (two wafls, say) can vary; one has realizations of grand canonical en- 
sembles. Obvious candidates to which the 0{n) model with symmetry breaking quadratic 
boundary terms studied in this work might apply are magnetic systems. Unlike fluids, 
magnetic systems do not lend themselves to direct or indirect types of measurements 
of thermodynamic forces of the kind used for fluid systems. Nevertheless, finite-size and 
residual free energies (as well as their temperature derivatives) are observable quantities 
that in principle can be measured. As far as Monte Carlo calculations are concerned, we 
do not see qualitatively new challenges. The simulation techniques that are currently in 
use for stTidying lattice n-vector models shoiild also work when easy-axes spin anisotropics 
in the boundary layers are included. Of course, in cases where some surface bonds are 
critically enhanced, careful — and cumbersome determinations of the corresponding 
critical enhancements are required. Wc hope that the present work will stimulate such 
Monte Carlo investigations, as well as experiments. 
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Appendix A. Calculation of required integrals 

In this Appendix we compute the integrals /}^''^''(L; f = 0) and /^^"'^"'^"'^"^(/.j f = 
0) defined by Eqs. (4.16) and (4.12), respectively. 

To compute i[^'^\L;0), note that the spectral decomposition of the free propagator 
GL,4>,a{P', z, z) is given by Eq. (3.6), except that we must leave out the r = sum- 
mand with eigenvalue fcg = when (AcB^) = (N,N). Substituting this expression into 
Eq. (4.16), we can perform the integration over p and exploit the fact that the eigen- 
functions ¥'r'^'^^(C) given in Eq. (3.7) are orthonormal on the interval [0, 1] to obtain 

, Ol-d„(3-d)/2 r3-d , 

where the prime on J^'r indicates omission of the r = summand when (Aa,Ba) = 
(N,N). 

Since the eigenvalues depend on the boundary conditions (A^, Bq), so do the series 
One finds 



1 for(A„,B„) = (N,N), 

1 for(A„,B„) = (D,D), (A.2) 

2^-'' - 1 for (A„, B„) = (D, N) or (N, D) , 



which together with Eq. (A.l) yields the results stated in Eq. (4.17). 
An analogous calculation yields 
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where we used the fact that 



Jo 

1+2 ^^1^2 for a = /3, 



= < 1 for (A^, ^ (A„, B„) = (D, N), (N, D), (A.4) 

1^2 ^^1^2 ^""^ Ac = Bq; = N and A/j = B/3 = D. 

The summations and can be carried out in a straightforward fashion. As 
results one obtains 

,(A,A;B,B) L^-2d2l-2d^d-Z [3^2(3 _ ^) ^ (3^^^ _ 1)^(5 _ 2rf)] 

^2 r2[(d-l)/2]cos2(dV2) ' ^ ^ ^ 



,(D,N;D,N)., ^ x^-2d2^-^-'7r''-3 [2(23-^^ - 1)^(^(3 - d) + {2^^^^ ^ 1)C(6 - 2d)] 
^ ^ ' ^ r2[(d-l)/2]cos2(d7r/2) 

(A.6) 



and 



(^'0)-J2 (L,0)- r2[(rf-l)/2]cos2(dV2) ■ ^ ^ ^ 

We need the e expansions of these integrals only to zeroth order. They read 

4^'^-'^'^\L;0) = ^^^+O{e), A,B = N,D, (A.8) 



and 



7f'™(L;0) = ^ + O(e), (A.9) 
4-^--'-'-\L-,0) = ir-'-'-\L-0) = +0(e). (A.10) 



Appendix B. Derivation of the representation (3.9) of g'"^'^\xi;x2) 

Here we give a derivation of Eq. (3.9) based on Poisson's summation formula (3.10). 

Wc start from the spectral representation (3.6), use this summation formula, and shift 
the integration variable k by 7r/2i. We thus arrive at 

The k integration can be performed. This gives 
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(-XJ „ 

4—-^ Jp 



(d-1) 



exp(ip ■ yia) \_\2Lj+z^-z,\./^ 
2 + ^ L 



g-|2Lj+02+zi|\/p^+f 



(B.2) 



When the result is expressed in terms of the free massive bulk propagator 



Gt\x\T)= I 
Jp 




(B.3) 



the representation (3.9) of ^ '{xi;x2) follows. 
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